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John Baker' ), Giovanni Vignale^ and A. G. Rojo^) 
^Department of Physics, University of Michigan, Ann Arbor, MI 48109-1120 
^Department of Physics, University of Missouri, Columbia, MO 65211 

In this paper we consider the effect of disorder on the non-dissipative Coulomb drag 
between two mesoscopic metal rings at zero temperature. Ring I has an Aharonov-Bohm 
flux present which creates a persistent current Jo. Ring 2 interacts with ring 1 via the 
Coulomb potential and a drag current, Jd is produced. We show that this drag current 
persists with finite disorder in each ring, and that for small disorder, Jd decreases with the 
square of the disorder amplitude. 

I. INTRODUCTION 

Electron-electron (e-e) interactions are responsible for a multitude of fascinating effects in condensed 
matter. They play a leading role in phenomena ranging from high temperature superconductivity and the 
fractional quantum Hall effect, to Wigner crystalization, the Mott transition and Coulomb gaps in disordered 
systems. The effects of this interaction on transport properties, however, are difficult to measure. A new 
technique has recently proven effective in measuring the scattering rates due to the Coulomb interaction 
directly^. 

This technique is based on an earlier proposal by Pogrebinsknu'u. The prediction was that for two con- 
ducting systems separated by an insulator (a semiconductor-insulator-semiconductor layer structure in par- 
ticular) there will be a drag of carriers in one film due to the direct Coulomb interaction with the carriers 
in the other film. If layer 2 is an "open circuit" , and a current starts flowing in layer 1 , there will be a 
momentum transfer to layer 2 that will start sweeping carriers to one end of the sample, and inducing a 
charge imbalance across the film. The charge will continue to accumulate until the force of the resulting 
electric field balances the frictional force of the interlayer scattering. In the stationary state there will be an 
induced, or drag voltage Vd in layer 2. 

There is a fundamental difference between transresistance and ordinary resistance insofar as the role of 
the Coulomb interaction is concerned. For a perfectly pure, translationally invariant system, the Coulomb 
interaction cannot give rise to resistance since the total current commutes with the Hamiltonian H . This 
means that states with a finite current are stationary states of H and will never decay, since the e-e interaction 
conserves not only the total momentum but also the total current. (For electrons moving in a periodic 
lattice, momentum and velocity are no longer proportional and the current could in principle decay by the 
e-e interaction.) If the layers are coupled by the Coulomb interaction, the stationary states correspond to a 
linear superposition of states in which the current is shared in different amounts between layers: the total 
current within a given layer is not conserved and can relax via the inter-layer interaction. n 

This mechanism of current degrading was studied in the pioneering experiment of Gramila et aln for GaAs 
layers embedded in AlGaAs heterostructures. The separation between the layers was in the range 200-500A. 
The coupling of electrons .and holes and the coupling between a two dimensional and a three dimensional 
system was also examinedcl. 

If we call I the current circulating in layer 1, the drag resistance (or transresistance) is defined as 

V D 

Pd = —- 

Most of the experiments done so far indicate the vanishing of pr> at zero temperature, something expected 
in the usual scattering theory of transport. _ 

The possibility of a drag-effect at zero temperature was considered by Rojo and MahanE3, who considered 
two coupled mesoscopidla rings that can individually sustain persistent currents, see Figure ([!]). The 
mechanism giving rise to drag in a non-dissipative system is also based on the inter-ring or inter-layer 
Coulomb interaction, the difference with the dissipative case being the coupling between real or virtual 
interactions. One geometry in which this effect comes to life is two collinear rings of perimeter L, with a 
Bohm-Aharonov flux, $i, threading only one of the rings (which we will call ring one). This is of course a 
difficult geometry to attain experimentally, but-has the advantage of making the analysis more transparent. 
Two coplanar rings also show the same effectEJ. If the rings are uncoupled in the sense that the Coulomb 
interaction is zero between electrons in different rings, and the electrons are non-interacting within the 
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rings, a persistent current Jq — —cdE/d^i = evp/L will circulate in ring onc0. If the Coulomb interaction 
between rings is turned on, the Coulomb interaction induces coherent charge fluctuations between the rings, 
and the net effect is that ring two acquires a finite persistent current. The magnitude of the persistent 
drag current Jd can be computed by treating the modification of the ground state energy in second order 

(2) 

perturbation theory AE , and evaluating 



Jd 



dAE, 



(2) 



d$ 2 



(1) 



$2=0 



with $2 an auxiliary flux treading ring two that we remove after computing the above derivative. 




FIG. 1. Schematic depiction of the non-dissipative drag setup. A persisent current Ji is induced in ring 1 by a 
Bohm-Aharonov flux. The Coulomb interaction V couples the charge fluctuations and generates a current in the 
second ring 

The question of the effect of disorder on persistent currents remains controversial. Since our project 
involves calculating the effect of disorder on an induced persistent current, we expect our results to shed 
some light on this issue. For an isolated pure ring the persistent current is Jo = cvf/L with L the perimeter 
of the ring and Vf the Fermi velocity. The most immediate effect of disorder is to introduce a mean free 
path t. One expects disorder to decrease the persistent current, and qualitative arguments indicate that it 
is decreased by a factor t/L: Jq — > evp/L(i/L). Our results indicate on firmer theoretical grounds that a 
similar argument can be used for the drag persistent current. 

In this paper we outline our detailed studies of the effect of disorder on non-dissipative drag using both 
analytic and numerical methods. 



II. GENERAL REMARKS ON THE NON-DISSIPATIVE DRAG 



The zero drag current can be finite only if quantum coherence, or entanglement, between the wave functions 
of the two systems is established. In this situation, the meaningful description of the dynamics of the 
combined system involves a single wave function, which distinghuishes from ordinary dissipative drag, a case 
in which one has scattering between two incoherently coupled systems. Figure ([!]) is a schematic illustration 
of this coherent coupling mechanism. 
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FIG. 2. Schematic depiction of the non-dissipative drag mechanism. A Bohm-Aharonov flux $ applied in ring 
1 displaces the k values of the free one-particle states by 4>/L, with 4>2ir$ / 4> , 4>o being the flux quantum. The 
unperturbed energy dispersions are therefore asymmetric in ring 1 and symmetric in ring 2. When the interaction is 
turned on, it creates virtual excitations of momentum — q in ring 1 and momentum q in ring 2 (shown by the dashed 
arrows). The amplitude of this excitation is V(q) /(Ai + A2). These excitations are not cancelled by those of reversed 
momentum (shown by the continuous arrows), the amplitude of which is V(q)/(Ai + A3) because of the asymmetry 
in the spectrum of 1. As a result there is a persistent current in ring 2 proportional to \V(q)\ 2 (f>/L from the terms 
indicated in this figure. 

We consider first two one-dimensional systems. Assume that, in the absence of the Coulomb coupling, 
system 1 carries a finite equilibrium current, which could in principle be established by an Aharonov-Bohm 
flux threading system 1 only. If system 2 is a one dimensional wire of perimeter 2irL, the mesoscopic nature 
of the zero drag current can be proven by the following analysis. 

Let ^>q be the ground state of the combined system. This wave function involves the coordinates of both 
systems. Let us consider system 2 as a closed ring geometry, and designate the coordinates of the particles 
in this subsystem as angular variables 9i, with i = 1, • • • , N 2 , and N 2 being the number of particles at system 
2. The kinetic component of the Hamiltonian of system 2 can then be written as 

H K - 2mL 2 2^ Q02 ■ W 

Consider the modified wave function constructed by applying a "boost" , or gauge transformation, on 



3 



the coordinates of system 2: 



N 2 

& = U(a)* = exp(ia^6 i )* , (3) 

i=l 

with a a parameter. By the variational theorem E' = (^>'\H\^>') > Eq, with H the Hamiltonian of the 
combined system, and Eq the total energy. On the other hand, explicit evaluation of E' gives 



with the current operator for system 2 given by 



JV 2 



27rmL 2 ^ 89. 



1=1 



Due to the variational nature of the bound, the dragged current has to obey the inequality: 

Jd rag M*ol4>o)<a 2 2^, (6) 

with p the particle density. Equation (^) emphasizes the mesoscopic nature of the dragged current: in the 
limit of L — » oo, Jdrag —> with the same length dependence as the persistent current in mesoscopic rings, 
the value of which is evF / L in the ballistic regime. Note that the bound is valid for strictly one-dimensional 
systems. 

Having established a bound, one needs to show that there is indeed a finite dragged current, and provide 
a quantitative estimate. We first present such a calculation treating the Coulomb interaction between the 
systems in second order perturbation theory. Consider two identical one-dimensional wires. Wire 1 is 
threaded by a Aharonov-Bohm flux 0i (in units of the flux quantum). In order to evaluate the induced 
current J2, we impose also a flux 02 in system 2, and compute 



j _ e dE 



(7) 

2=0 



We neglect the Coulomb interaction within each wire, and consider the ballistic regime (no impurities in 
either system). In the absence of coupling, and for both fluxes 0$ < ir/2, the ground state consists of two 

Fermi systems with one particle energies E^ — ^L L i (jhi — (f>i) 2 , and occupied levels for rii < rip (i = 1, 2, 
and rip = N/2, N being the particle number at each ring). Let V(q) be the Fourier transform of the 
Coulomb coupling, which for wires separated a distance d has the form V(q) = (2e 2 /L)K (qd), K (x) being 
the zero-order Bessel function of imaginary argument. The second order correction to the energy is then 
given by: 

= _mTJ v V 2 (g)/n 1 (l-/n 1+ q)/n 2 (l-/n 2 -q) (R] 
2 h 2 n ^ Q (Q + ni + 0i - n 2 - 0a) ' { 1 

,712 

with Q,ni,ri2 integers, and f m Fermi functions: f m = 1 if |m| < rip, and zero otherwise. The above sum 
is now evaluated transforming the sum into integrals over the continuum variables q = Q/L, ki = rii/L. 
Evaluating the integrals, and computing the derivative with respect to 02, we obtain 

me 5 1 I(kpd) 

with I(kpd) = dq q ^ . In the limit of large kpd, which corresponds to the intcrparticle distance being 
much smaller than the distance between the systems, we obtain 

J2 " J °(Jao) 2 (Jd) 2 ' (10) 
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with Jo = —evpcfri/ L being the persistent current carried by the otherwise uncoupled system 1, and ciq the 
Bohr radius. We have proven that there is an induced persistent current due to the Coulomb interaction. 
We now ask ourselves about the induced effect if system 2 is made open, so that no current can circulate. 
In the transport situation, a voltage will be induced. Here, we show that there is no voltage induced. We 
start with a setup that, in the absence of the flux in system 1, is "parity even". By this we mean that the 
charge distribution in wire 2 is symmetric around the center of the wire. We want to know if this symmetry 
is broken by applying the flux in system 1, an operation that breaks the time reversal symmetry. Let us call 
P and T the parity and time reversal operators that interchange the ends of the wire. We want for example 
the induced dipole moment in wire 2, x 2 = (^ol^l^o)- The operator PTx2(PT)~ 1 — —X2, while the wave 
function is invariant under PT, which implies Xi = 0, hence there is no induced voltage. 



III. DISORDER AND THE NON-DISSIPATIVE DRAG 



In this section we outline our results on the effect of disorder in the non-dissipative drag. In calculating 
the effects of disorder we use the two ring geometry considered by Rojo and Mahan, see Figure (|l|), and 
calculate the second order Coulomb interaction between the conduction electrons in the two rings. The 
Coulomb potential is 

V = Y2v k pk,iP-k,2 = J dx J dx'p 1 (x)p 2 (x')V(x - x'), 

With pi the charge density at ring i. The second order correction to the ground state energy due to the 
Coulomb interaction is 

A£ _ V V I E fc VkU*Pne™ k */ L r n/ dx)(f xj> m j* k *'l L r m .dx')f n >{l - UMl - f m )\ 2 , , 

^ ^ E n - E n , +E m - E^. 



n,n' ra^vn' 



where ip n is an eigenstate in the presence of disorder. From this expression for the energy shift we can 
calculate the drag current from Equation (flh. 



A. Analytics 

In this section we estimate the effect of disorder on the non-dissipative drag current for the case in which 
disorder is present only in the ring on which the Bohm-Aharonov flux is applied. The driven ring (ring 2), 
on which the drag current circulates, will be taken as disorder-free. Momentum remains a good quantum 
number in ring 2 making the calculation more tractable. The first order correction to the wave function is 
given by 

q k v H ' 

where E^ are the one-particle energies for the states of ring 2, and \v) is a many-body state of ring 1 with 

energy E e . The ground state of ring 1 is I^Jq 1 ^), and its energy is £?o,i- Now, since we are neglecting 
interactions within each ring, the resulting equilibrium current in ring 2 is given by 

t e ^hq /fc (l-/ fc+9 )/ p (l-/,)|(^>)| 2 

J2 ^V- rh ^ q -E k+Ev -E^ ' (13) 

where now \v) refers to the exact one-particle states with energies E v corresponding to the disordered 
Hamiltoninan in ring 1. We can rewrite the above espression in terms of the spectral function S(q,Lo) 
defined as 

S(q,u) = - U)\(p\e l **\v)\ 2 6 [w - {E v - E u )/h) . (14) 
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We will consider the function S(q,uj) in the approximation in which the matrix element is given 

by the diffusive lorentziantLil: 

M e "»| 2 = nhN(0) WY + {E,-E v y/h^ (15) 

where D is the difussion constant and N(0) is the density of states of the system. In this approximation we 
obtain that S(q, u>) is given by 

S( q ,u) = N(0) {D ^/ +u;2 . (16) 

Before replacing this expression in Equation ([l3]) let us recall that there is-a flux <E> threading ring 1 and 
therefore one expects S(q,tj) S(—q,u>). We follow Ambegaokar and EckerrJllJ in including the effect of the 
flux in the diffusive motion through the replacement: 



Dq*^D<? = Dlq-wj- ) , (17) 

with <j) being the flux in units of the flux quantum. 
The induced current will therefore be given by 

J 2 = ty^^-\V(q)\*hDfN(0)Y: [ du .J^-J^l .„ " .. . (18) 
£Y m k' {E k+q - E k + hLu) 2 {Dq 2 +uj 2 ) 



For small wavevectors {q <C fcf ) we have: 



\ - /fc(l ~ fk+g) = L q 



k 



and also, in the limit of ql < 1, with I being the mean free path 

1 uj (ql) 

(Dq 2 ) 2 +uj 2 ~ (hv F ) 2 q 



duj J~ 2 ^2(n*2V i -.2 - 71—^272- ( 2 °) 

(^9 + ^1 



We are interested in the lowest order in <fi for the induced current, which gives 



i-i^jEW. < 21 » 



47r mu4 L 



which we can now rewrite using D = v F t as 

£ N(0){e 2 /d) 2 C 



J: 

where C is a constant, 



x 0, (22) 



C = / dxx 2 K {x) 2 = .308425 (23) 
Jo 

The first term in square brackets in Equation ( |22| ) corresponds to a familiar expression for the persistent 
current in ring 1 in the presence of disorder. The value of terms in the second square bracket can be computed 
taking N(0) = 1/A, with A ~ lOi^T being the level spacing for a ring of L ~ lfim, Ep — 2eV, and a distance 
between rings of d = 100A. Note that this term contains the product of two ratios: a small one given by 
Eq ou \/E f , with Seoul = e2 /d, and a large one given by £coui/A. This gives a number of order one, a result 
that probably overestimates the drag current, but serves as an indication that the effects of disorder are not 
extreme. The second square bracket also contains an additional ratio, the mean free path to the distance 
between rings. This additional factor shows that the effects of disorder are stronger in the drag current from 
that in the driving ring. In order to test this results we performed numerical simulations, which we present 
in the following sections. 



(i 



B. Numerical simulations 



1. Perturbative treatment of the Coulomb interation 



In evaluating the drag current computationally we consider a discrete ring with N lattice sites and P < N 
electrons. We model disorder by placing a random disorder potential at each lattice site. The hamiltonian 
for an electron hopping between lattice sites in this ring is given by 

N N 

H = t(J2 CjC t+1 e^ + Ct-xCie-*) + £ W n ClC n , 

i—l i n—1 

where <j> is the magnetic flux through the ring, C\ is the electron creation operator at site i and w n is the 
disorder potential at site n. For N lattice sites, this gives an N x N hopping matrix. 

In computing the energy shift for the two ring system we work with the x space representation of Equa- 
tion©, 

AE _ y, y, I Ex Ex' V ( x - x ') < x \ n >< n '\ x >< x '\ m >< m '\ x ' > fn'jl - fn)fm'(l ~ fro)? 



, E n — E n / + E m — E n 

n,n — 1 m,m — 1 



(24) 



Here x and x' denote discrete positions of the lattice sites in rings one and two respectively and the \n >'s 
and -En's are the eigenvectors and eigenvalues obtained numerically from the hopping matrix. We obtain 
disorder averaging by evaluating AE with different realizations of the random disorder potentials, W n , at 
values between — W and W where W is the disorder amplitude. The result of the computer simulations are 
shown in figure @ for a system of 10 lattice sites and 7 particles. The ratio of the drag current to its zero 
disorder value Jd/Jd{0) is plotted both for a system in which disorder is present in ring 2 only and for a 
system of two disordered rings. The ratio J o /J o (0) is also plotted. For small disorder amplitude, Jd oc W 2 . 



7 



0.0 1.0 2.0 3.0 4.0 

w/t 



FIG. 3. Plot of Drag Current vs. Disorder Amplitude for a system of 10 lattice sites and 7 particles with disorder 
averaging. We show three curves: Jd/Jd(0) for one disordered ring (square symbols), Jd/Jd(0) for both rings with 
disorder (triangles) and J o /J o (0) (circles). 

2. Non-perturbative treatment for very small rings by Lanczos method 

In this section we present some exact results for small clusters. We use the Lanczos method to diagonalize 
the problem, and obtain results that are non-perturbative in the interaction. As a first illustration, Figure (^) 
shows the persistent and drag currents both with and without disorder. The drag current follows the 
persistent current of ring 1 in its periodicity of one flux quantum as a function of the applied flux through 
ring 1. 
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Persistent and drag current vs. AB flux 

(with and without disorder) 

i i i i i i i i i i i i i i i 





FIG. 4. Exact results for the persistent current and the drag current as a function of flux for two rings of six sites, 
each of them with two particles. We show curves for zero disorder as well as finite disorder. 

Figure (||) shows the drag current for two systems of different sizes. Note that the dependence with 
disorder is stronger for the larger system as expected from the factors of £/L that appear in the analytical 



expressions in section III A 
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FIG. 5. Exact results for the drag current as a function of disorder amplitude W in units of the hopping matrix 
element for two rings interacting via a delta function potential of amplitude 0.5i. The flux in ring 1 is $1 = O.20o 

In conclusion we have established that the drag current remains finite for finite disorder. We have shown 
by numerical simulations of finite clusters and by analytical considerations that the effect of disorder on the 
drag current is more pronounced than the effect of disorder on the persistent current in a single ring. 
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